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We present accurate predictions of the effective Majorana mass \mpp\ in neutrinoless double-/? 
decay in the standard case of 3 v mixing and in the case of 3+1 neutrino mixing indicated by the 
reactor, Gallium and LSND anomalies. We have taken into account the uncertainties of the neutrino 
mixing parameters determined by oscillation experiments. It is shown that the predictions for \mpp\ 
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the literature, and that future measurements of neutrinoless double-/? decay and of the effective 
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I. INTRODUCTION 

Neutrino flavor oscillations have been observed in 
many solar, reactor and accelerator experiments (see the 
recent reviews in Refs. BED, in agreement with the cur¬ 
rently standard paradigm of three-neutrino (3+ mixing. 
The global fits of neutrino oscillation data in the frame¬ 
work of 3v mixing [3H5] give us rather precise information 
on the values of the elements of the three mixing angles 
which parameterize the neutrino mixing matrix and on 
the values of the two independent neutrino squared-mass 
differences, the smaller “solar” squared-mass difference 
Am| 0 L ~ 7.5 x 10 _5 eV 2 and the larger “atmospheric” 
squared-mass difference Am^ TM ~ 2.4 x 10 _3 eV 2 . 

However, the standard 3v mixing paradigm has been 
challenged by indications in favor of short-baseline os¬ 
cillations generated by a new larger squared-mass differ¬ 
ence Am| BL ~ leV 2 : the reactor antineutrino anomaly 
[B] , which is a deficit of the rate of v e observed in several 
short-baseline reactor neutrino experiments in compar¬ 
ison with that expected from the latest calculation of 
the reactor neutrino fluxes E0; the Gallium neutrino 
anomaly [9ifl3]. consisting in a short-baseline disappear¬ 
ance of v e measured in the Gallium radioactive source 
experiments GALLEX [14] and SAGE [15]; the signal 
of short-baseline —» v e oscillations observed in the 
LSND experiment [16] [17] . The simplest extension of 
3v mixing which can describe these short-baseline oscil¬ 
lations taking into account other constraints is the 3+1 
mixing scheme ii, in which there is an additional 
massive neutrino at the eV scale and the masses of the 
three standard neutrinos are much smaller. Since from 
the LEP measurement of the invisible width of the Z bo¬ 
son we know that there are only three active neutrinos 
(see Ref. [20]), in the flavor basis the additional massive 
neutrino corresponds to a sterile neutrino [21] , which does 
not have standard weak interactions. 


A fundamental questions that remains open is: are 
neutrinos Dirac or Majorana particles? This question 
cannot be investigated in neutrino oscillation experi¬ 
ments, where the total lepton number is conserved and 
there is no difference between Dirac neutrinos with a con¬ 
served total lepton number and truly neutral Majorana 
neutrinos, which do not have a conserved total lepton 
number. The most promising process which can reveal 
the Majorana nature of neutrinos is neutrinoless double¬ 
beta decay, in which the total lepton number changes by 
two units (see the recent review in Ref. [22]). 

From the present knowledge of the neutrino squared- 
mass differences and mixing angles it is possible to pre¬ 
dict the possible range of values for the effective Ma¬ 
jorana mass \mpp\ in neutrinoless double-beta decay as 
a function of the absolute scale of neutrino masses (see 
Ref. [22]), which is still unknown, up to the upper bound 
of about 2 eV at 95% C.L. established by the Mainz [23] 
and Troitsk [21] Tritium /3-decay experiments. 

The introduction of a sterile neutrino at the eV mass 
scale can change dramatically the prediction for the pos¬ 
sible range of values for the effective Majorana mass in 
neutrinoless double-beta decay [13| 1251133] . In this paper 
we present accurate predictions for \mpp\ taking into ac¬ 
count the results of the global fit of solar, atmospheric 
and long-baseline reactor and accelerator neutrino oscil¬ 
lation data presented in Ref. [3] and the results of an 
update [34] [35] of the global fit of short-baseline neutrino 
oscillation data presented in Ref. [TJJ]. We are partic¬ 
ularly interested to determine accurately the conditions 
for which \mpp\ > 0.01 eV, which may be probed exper¬ 
imentally in the near future (see Refs. [2BMT] ). and the 
conditions for which there can be a cancellation of the 
different mass contributions to \mpp\, which leads to an 
unfortunate uncertainty for the possibility of ever observ¬ 
ing neutrinoless double-beta decay (unless it is induced 
by new interactions and/or the exchange of new particles; 
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see Refs. 

The plan of this paper is as follows. In Section [II] 
we discuss the predictions for \mpp\ in the standard 3v 
framework, taking into account the two possible normal 
and inverted mass orderings. In Section |TTT| we discuss 
how these predictions are modified in the 3+1 mixing 
framework. In Section HVl we draw our conclusions. 


II. THREE-NEUTRINO MIXING 

In the standard three-neutrino (3^) mixing framework, 
the effective Majorana mass in neutrinoless double-beta 
decay is given by 

\ m pp\ = |#i + M 2 e m2 + /i 3 e* a3 1 , (1) 

where 

= |£^efc| (2) 

is the partial contribution of the massive Majorana neu¬ 
trino Uk with mass mk- The elements U e k of the mix¬ 
ing matrix, which quantify the mixing of the electron 
neutrino with the three massive neutrinos, can have un¬ 
known complex phases, which generate the two complex 
phases a 2 and a 3 in Eq. 0 . Since the values of these 
phases is completely unknown, all predictions of the value 
of \mpp\ must take into account all the possible range of 
these phases, from 0 to 2n. 

We use the results of the global fit of solar, atmospheric 
and long-baseline reactor and accelerator neutrino oscil¬ 
lation data presented in Ref. [3] , which are given in terms 
of the mixing angles $ 12 , #13 that determine the abso¬ 
lute values of the first row of the mixing matrix U in the 


standard parameterization: 

\U e l\ = COS #13 COS #12 , (3) 

| Ue 2 1 = COS #13 sin #12, ( 4 ) 

\U e3 \ = sin#i 3 . ( 5 ) 

The results for the neutrino squared-mass differences are 
expressed in terms of the solar and atmospheric squared 
mass differences, which are defined by 

a ™sol = Ara2!, (6) 

Am ATM = 2 l Am 31 + Am 32| > (?) 


where = m^ — m|. Given this assignment of the 

squared mass differences, it is currently unknown if the 
ordering of the neutrino masses is normal (NO), such 
that mi < m 2 < m 3 or inverted ( 10 ), such that m 3 < 
mi < m 2 . We discuss these two cases separately in the 
following subsections. 


A. Normal Ordering 

In order to study the case of Normal Ordering (NO), 
we express the neutrino masses in terms of the lightest 



FIG. 1. Best-fit values (b.f.) and la, 2a and 3a allowed 
intervals of the three partial mass contributions to \mpp\ in 
Eq. 0 as functions of the lightest mass mi in the case of 3v 
mixing with Normal Ordering. 



FIG. 2. Ratio m^J/m^ (see Eq. (13)) as a function of mi 
for the best-fit values of the partial mass contributions in the 
case of 3u and 3+1 mixing with Normal Ordering of the three 
lightest neutrinos. 
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FIG. 3. Value of the effective Majorana mass \mpp\ as 
a function of the lightest neutrino mass mi in the case of 
3 v mixing with Normal Ordering. The signs in the legend 
indicate the signs of e 1012 , e tc * 3 = ±1 for the four possible cases 
in which CP is conserved. The intermediate yellow region is 
allowed only in the case of CP violation. The 90% upper limit 
is explained in the main text. 


neutrino mass m m i n : 


mi = 

^min? 

(8) 

m 2 = 

\/ m min + Am SOL> 

(9) 

m 3 = 

J m min + Am ATM + A "»SOl/ 2 - 

(10) 


Figure [I] shows the best-fit values and the lcr, 2a and 
3cr allowed intervals of the three partial mass contribu¬ 
tions to \mpp\ in Eq. (jTj) as functions of the lightest mass 
mi. We calculated the confidence intervals using the x 2 
function 


xiu — X 2 ( / ^ m soh) + X 2 (^ m ATM) 

+ X 2 (sin 2 tf 12 ) +x 2 (sin 2 tf 13 ), (11) 


TABLE I. Ranges of mi, mp and E for which there can be a 
complete cancellation of the three partial mass contributions 
to \mpp | for the best-fit values (b.f.) of the oscillation pa¬ 
rameters and at lcr, 2cr and 3cr in the case of 3v mixing with 
Normal Ordering. 

b.f. lcr 2 <j 3a 

mi[10“ 3 eV] 2.3-6.6 1.9 - 7.2 1.6 - 8.0 1.3 - 9.0 

mp [10 -2 eV] 0.9 - 1.1 0.9 - 1.2 0.8 - 1.3 0.8 - 1.4 

E [HT 2 eV] 6.1 - 6.8 5.9 - 7.0 5.8 - 7.2 5.7 - 7.4 


FIG. 4. Value of the effective Majorana mass \mpp\ as a 
function of effective electron neutrino mass mp in Eq. (22) in 
the case of 3u mixing with Normal Ordering. The legend is 
explained in the caption of Fig. [3] The limits are explained 
in the main text. 



E [eV] 


FIG. 5. Value of the effective Majorana mass \mpp 


a function of sum of the neutrino masses E in Eq. (23) in 


the case of 3v mixing with Normal Ordering. The legend is 
explained in the caption of Fig. [3] The limits are explained 
in the main text. 


with the partial x 2 ’ s extracted from Fig. 3 of Ref. [3], 
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neglecting possible small correlations of the four mixing 
parameter^] For each value of m\ we calculated the con¬ 
fidence intervals for one degree of freedom. Although this 
method is in principle better than the method based on 
the propagation of errors of the parameters used by many 
authors, in practice it leads to similar results, because 
the x 2, s of Amg OL , Am^ T M) sin 1 2 $12 and sin 2 $13 in 
Fig. 3 of Ref. [3] are very well approximated by quadratic 
functions, which correspond to Gaussian uncertainties for 
which the method of propagation of errors is valid. 

From Fig. [l] one can see that for mi < 2 x 10 -3 eV 
the contribution of fi 2 is dominant and cannot be can¬ 
celed by the smaller contributions of fi\ and fi 3 for any 
values of the relative phase differences. In the inter¬ 
val m 1 « (2 - 7) x 10 _3 eV cancellations are possible, 
mainly between fi\ and /12 which have similar values, with 
the smaller contribution of / 13 , which is about 2.3 times 
smaller than p- 2 - For mi > 7 x 10 _3 eV again there can¬ 
not be a complete cancellation, because the contribution 
of /ii is dominant. 

The result for \mpp | is shown in Fig. [3j where we have 
plotted separately the allowed bands for the four possible 
cases in which CP is conserved ( 0 : 2 , 0 L 3 = 0, 7 r) and the co¬ 
efficients of the contributions are real. These are the ex¬ 
treme cases which determine the minimum and maximum 
values of \mpp\. The areas between the CP-conserving 
curves correspond to values of \mpp\ which are allowed 
only in the case of CP violation [31] 1501154] (although 
there is no manifest CP violation |55J). 

Figure [3] shows also the 90% C.L. upper limit band 
for \m$$\ estimated in Ref. [22] from the results of the 
KamLAND-Zen experiment [56] taking into account the 
uncertainties of the nuclear matrix element calculations. 
The reliability of this upper limit is supported by the 
upper limits with the same order of magnitude following 
from the results of the Heidelberg-Moscow [57], IGEX 
[58], GERDA [59], NEMO-3 [60], CUORICINO [0], and 
EXO [HI] experiments. 

From Fig. [3] one can see that, in agreement with the 
discussion above, there can be a complete cancellation of 
the three partial mass contributions to \mpp\ for mi in 
the intervals given in Tab. [I] at different confidence levels. 

The exact determination of the region in which there 
can be a complete cancellation of the partial mass contri¬ 
butions to \mpp\ in the general case of IV-neutrino mixing 
can be done in the following relatively simple wa}0 For 
each value of the lightest mass mi let us denote by a the 
index of the largest mass contribution, i.e. 


Then, we can consider the quantities 

= n a ±^2fik- (13) 

k^a 

The quantity is always positive and represents the 
most favorable case, in which all the mass contribution 
add with the same phase, giving the maximum value of 
\mpp\ for any value of the unknown phases: 

|^/3/3|max ^/ 3/3 ‘ (14) 

The quantity m^ represents the extreme case in which 
the phases of all the other partial mass contributions are 
equal and opposite to the phase of the largest mass con¬ 
tribution. It is evident that if rrr^ > 0, the value of 
mpp gives the minimum possible value of \mpp\ for any 
value of the unknown phases, because it corresponds to 
the maximal cancellation between fi a and the maximum 
Yhkj^a Mfc °f the other partial mass contributions. On the 
other hand, if 

”47 < o, (15) 

there is an intermediate value of the phases which gives 
\mpp\ = 0. This can be seen clearly by writing \mpp\ as 


with 


l m /?/?l = 


Mo + e 


ia! 



(16) 


/' 


! 


k^a 


(17) 


with an unknown phase a' and N — 2 unknown phases £&, 
of which one can be fixed to zero. The only possibility 
to have \mpp\ = 0 can be realized for of == tt if y! = Mo- 
This equality can occur only if 


< /-'« < lAn 


(18) 


where Mmin an< l Mmax are th e minimum and maximum 
values of fi! for any value of the unknown phases Since 
we have always Mmin — Mo because of Eq. (12) and 


= 

k^a 


(19) 


the inequality in Eq. ( [15] ) is the necessary and sufficient 
condition for having \mpp\ = 0 for some value of the 
unknown phases. Hence, the minimum value of \mpp\ is 
given by 


Mo > Mfc for k^a. 


(12) 1771,901 min = max 

which can also be written as 


"47° 


( 20 ) 


1 The only significant correlations discussed in Ref. [3] are those 
which involve the mixing angle $ 23 , which is irrelevant for neu¬ 
trinoless double-/? decay, and the Dirac phase, whose effect in 
neutrinoless double-/? decay is masked by the two unknown Ma- 
jorana phases. 

2 Other ways are discussed in Refs. [3011631164J . 


l m /?/?l min = max[2/t fe - \rripp\ 

max? 0]. (21) 


3 Equation ( |2l| ) is the generalization to iV-neutrino mixing of 
that obtained with a different proof in Ref. | 63| in the case of 3z/ 
mixing. 
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Fig. [ 2 ] shows the value of the ratio as a 

function of m\ for the best-fit values of the partial mass 
contributions, which is sufficient for the determination 
of the interval of mi for which there can be a complete 
cancellation of the partial mass contributions. One can 
see that in the case of 3 v mixing mffcp is negative and 
it is possible that \mpp\ = 0 only in the interval of mi 
given in Tab. [TJ 

Let us now consider the opposite possibility that \mp$\ 
is larger than about 0.01 eV, which is a value that may be 
explored experimentally in the near future. From Fig. [3] 
one can see that \mpp\ > 0.01 eV can be realized only 
for mi > 0.008 eV. This range of mi corresponds to 
almost degenerate m\ and m 2 , because ^/Am| OL « 8.7 x 
10 -3 eV. Hence, it will be very difficult to measure \mpp\ 
if there is a normal hierarchy of neutrino masses (mi <C 
m 2 ^ m 3 ) for any value of the unknown phases 012 and 

a 3 in Eq. ([TJ. 

One can also see from Fig. [3] that \mpp\ > 0.01 eV 
is realized independently of the values of the unknown 
phases 0 L 2 and ct 3 for mi > 0.04 eV, which is close to the 
region mi > yj AmJ TM « 0.05 eV in which all the three 
neutrino masses are quasidegenerate. 

Figure [3] gives a clear view of the possible values of 
\mpp\ depending on the scale of the lightest mass mi, but 
it is of little practical usefulness, because it will be very 
difficult to measure directly the value of mi. In practice, 
the investigation of the absolute values of neutrino masses 
is performed through the measurements of the effective 
electron neutrino mass 

mp = ^/|t/ e i| 2 m? + |( 7 e 2 | 2 ml + |C/ e 3 | 2 m| ( 22 ) 

in /3-decay experiments [23j [24] and through the mea¬ 
surement of the sum of the neutrino masses 


E = mi + m 2 + m 3 


(23) 


in cosmological experiments (see, for example, Ref. [65] 5. 
Hence, it is useful to calculate the allowed regions in the 
mp~\mpp\ and T,-\mpp\ planes [664469] . which are shown 
in Figs. [4] and [5] In this case, the confidence intervals are 
calculated using the y 2 function in Eq. (11) with two de¬ 


grees of freedom. We have plotted separately the allowed 
bands for the four possible cases in which CP is conserved 
(a 2 , a ?3 = 0, 7 r), in order to show the regions in which CP 
is violated. Potentially the possibility of measuring val¬ 
ues of \mpp\ and mp or E in these regions is very exciting 
for the discovery of CP violation generated by the Majo- 
rana phase^J but in practice such measurement is very 
difficult because it would require a precision which seems 


4 The phases a .2 and 0:3 in Eq. (MJ) depend on the values of one 
Dirac phase and two Majorana phases in the neutrino mixing ma¬ 
trix (see Ref. [70]). The Dirac phase can be measured in neutrino 
oscillation experiments and there is some indication on its value 
[3}f5]. On the other hand, the values of the two Majorana phases 



FIG. 6. Best-fit values (b.f.) and lcr, 2cr and 3cr allowed 
intervals of the three partial mass contributions to \mpp\ in 
Eq. 0 as functions of the lightest mass m 3 in the case of 3u 
mixing with Inverted Ordering. 


to be beyond what can be currently envisioned [7TU74] . 
especially taking into account the current uncertainty of 
the calculation of the nuclear matrix element in neutri¬ 
noless double-/? decay (see Refs. EatEHZZI). 

Figures [4] and [ 5 ] show the same 90% C.L. upper limit 
band for \mpp\ as in Fig. [ 3 J In addition, Fig.[i]show^]the 
90% C.L. sensitivity on mp of the KATRIN experiment 
[75] . which is scheduled to start data taking in 2016, and 
Fig.[5]shows the 95% bayesian upper limit on E obtained 
by the Planck collaboration [65] . 

The intervals of mp and E for which there can be a 
complete cancellation of the three partial mass contribu¬ 
tions to \mpp\ are given in Tab. [i] On the other hand, 
from Figs. [I] and [ 5 ] one can see that \mpp\ > 0.01 eV for 
any value of the unknown phases ct 2 and ck 3 for mp > 
0.05 eV and E > 0.15 eV. The 3cr lower bounds for mp 
and E are, respectively, 0.8 x 10 -2 eV and 5.6 x 10 -2 eV. 


can be measured only in lepton-number violating processes such 
as neutrinoless double-/! decay. In the future, if the Dirac phase 
will be measured, CP violation in neutrinoless double-/! decay 
may provide information on the Majorana phases. 

5 The most stringent current upper limits on mp obtained in the 
Mainz (mp < 2.3 eV at 95% C.L.) | 23| and Troitsk (mp < 2.1 eV 
95% C.L.) [24J experiments are out of the scale in Fig. [ 4 ] 
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m p [eV] 


FIG. 7. Value of the effective Majorana mass \mpp\ as a func¬ 
tion of lightest neutrino mass in the three neutrino case for 
the Inverted Ordering. The legend is explained in the c aptio n 
of Fig. [ 3 ] The 90% upper limit is explained in Section eai 


FIG. 9. Value of the effective Majorana mass Ira^l as a 
function of effective electron neutrino mass mp in Eq. (22) in 
the case of 3 v mixing with Inverted Ordering. The legend is 
explained in the caption of Fig. [3] The limits are explained 
in in Section HU] 
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FIG. 8. Ratio ( see Eq. (13)) as a function of 

m 3 for the best-fit values of the partial mass contributions in 
the case of 3u and 3+1 mixing with Inverted Ordering of the 
three lightest neutrinos. 
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FIG. 10. Value of the effective Majorana mass |mpp 


a function of sum of the neutrino masses E in Eq. (23) in 


the case of 3v mixing with Inverted Ordering. The legend is 
explained in the caption of Fig. [3] The limits are explained 
in in Section fll Al 
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B. Inverted Ordering 

In the case of Inverted Ordering (10), the expressions 
of the neutrino masses in terms of the lightest neutrino 
mass m m i n are: 


phase 0:4 that must be varied from 0 to 27 t as a 2 and a?, 
in order to calculate the predictions of the value of \mpp\. 

The absolute values of the relevant first row of the 4x4 
mixing matrix U is given by the simple extension of the 
standard parameterization: 


m 3 = 

Ulmin? 

(24) 

mi = 

\J m min + Am ATM “ Am SOL/ 2 > 

(25) 

m 2 = 

J m min + Am ATM + A ”*SOl/ 2 - 

(26) 


| U gi | = COS #14 COS #13 COS #12, (30) 

| U e 2 1 = COS #14 COS #13 sin #12, (31) 

| u e3 1 = cos #14 sin #13 , (32) 

\U e 4 \ = sin #14. (33) 


Figure [ 6 ] shows the best-fit values and the lcr, 2 cr and 
3cr allowed intervals of the three partial mass contribu¬ 
tions to \mpp\ in Eq. (JT|) as functions of the lightest mass 
m 3 . One can see that fi\ is always dominant, because 
#12 is smaller than 7 r /4 and \U e \\ > \U e 2 \ > \U e $\. There¬ 
fore, in the case of Inverted Ordering there cannot be a 
complete cancellation of the three mass contributions to 
\mpp\ (Fig. 8 shows that mi"^ is always positive) and we 
obtain frornTig. [7] the lower bounds 


Since in the case of 3+1 neutrino mixing, as well as in 
any extension of the standard Zv mixing, the ordering of 
the three standard massive neutrinos is not known, in the 
following two subsections we consider separately the two 
cases of Normal and Inverted Ordering of zq, zq, ^ 3 . The 
values of their masses as functions of the lightest mass 
ra m in are given by Eqs. (j 8 |-( 10) in the Normal Ordering 
and by Eqs. (24)-(26) in the Inverted Ordering. However, 
in both cases we have 


\mpp\ > 1.6 (lcr), 1.5 (2a), 1.3 (3a) x 10 “ 2 eV. (27) 

In the case of an Inverted Hierarchy (m 3 «mi < m 2 ) 
we also have the upper bounds 

\mpp\ < 4.8 (la), 4.9 (2a), 4.9 (3a) x 10“ 2 eV. (28) 


m 4 - # m rnin + Am §BL> ( 34 ) 

neglecting the contributions of Amg OL and Am\ TM , 
which are much smaller than Am| BL . We calculated the 
confidence intervals using the y 2 function 


The next generations of neutrinoless double-beta decay 
experiments (see Refs. [36h4T] ) will try to explore the 
range of \mpp\ between the limits in Eqs. (27) and (28), 
testing the Majorana nature of neutrinos in the case of 
an Inverted Hierarchy. 

Figures [9] and 10 show the correlation between \mpp\ 
and the measurable quantities mp and E in the Inverted 
Ordering. Since in this case both mp and E have rela¬ 
tively large lower bounds (4.6xlO _2 eV and 9.4 xlO -2 eV, 
respectively, at 3a) there is a concrete possibility that 
near-future experiments will determine an allowed region 
in these plots if in nature there are only three neutrinos 
with Inverted Ordering. 


III. 3+1 MIXING 

In this section we consider the case of 3+1 mixing 
in which there is a new massive neutrino z /4 at the eV 
scale which is mainly sterile. As explained in Section [TJ 
3+1 mixing is motivated mm by the explanation of 
the reactor, Gallium and LSND anomalies, which re¬ 
quires the existence of a new squared-mass difference 
AmgBL ~ leV 2 . In this case, the effective Majorana 
mass in neutrinoless double-beta decay is given by 

I'ffl'Pp | = |/ii + ^2^ ia2 + + A4e iQ:4 1 , (29) 

with the partial mass contributions given by Eq. The 
contribution of v 4 enters with a totally unknown new 


Xl+l = XL + X 2 ( Am SBL+ in2 ^14), 


(35) 


with xi v defined in Eq. (11) and x 2 ( Am g BL , sin 2 $ 14 ) 
obtained from an update [34N35] of the global fit of short- 
baseline neutrino oscillation data presented in Ref. [T9] . 

After Eq. ( pTj ) we noted that in the case of Zv mix¬ 
ing our statistical method for the calculation of the un¬ 
certainty of \mpp\ and the usual method based on the 
propagation of errors lead to similar results, because 
the x 2 ’ s °f th e relevant 3z/ mixing parameters are very 
well approximated by quadratic functions. On the other 
hand, the usual propagation of errors is inaccurate in 
the case of 3+1 mixing, because the marginal % 2 ’ s of 
AmgBL an d sin 2 #14 are not quadratic. Moreover, there 
are significant correlations between Am| BL and sin 2 #14 
(see Fig. 3 of Ref. [19]) which are taken into account in 
X 2 (AmgBL? si n 2 #14). 


TABLE II. Ranges of mi, mp and E for which there can be 
a complete cancellation of the four partial mass contributions 
to \mpp\ for the best-fit values (b.f.) of the oscillation param¬ 
eters and at la, 2a and 3a in the case of 3+1 mixing with 
Normal Ordering of the three lightest neutrinos. 

b.f. la 2a 3a 

mi[10“ 2 eV] 3.5 - 10.5 2.8- 12.5 2.0 - 16.3 1.3- 20.0 

mp [10 -2 eV] 3.6- 10.5 2.6- 14.3 1.9 - 17.8 1.5 - 22.7 

E [10~ 2 eV] 13.3 - 32.6 10.4 - 43.8 8.7 - 53.9 7.6 - 68.5 


















Lightest mass: rrq [eV] 



1CT 2 1CT 1 1 

m p [eV] 


FIG. 11. Best-fit values (b.f.) and la, 2a and 3a allowed 
intervals of the four partial mass contributions to \mpp\ in 
Eq. (29) as functions of the lightest mass mi in the case of 3+1 


mixing with Normal Ordering of the three lightest neutrinos. 


FIG. 13. Value of the effective Majorana mass \mpp \ as a 
function of effective electron neutrino mass mp in Eq. \22\ 
in the case of 3+1 mixing with Normal Ordering of the three 
lightest neutrinos. The legend is explained in the caption of 
Fig. |T2| 



FIG. 12. Value of the effective Majorana mass \mpp\ as a 
function of the lightest neutrino mass mi in the case of 3+1 
mixing with Normal Ordering of the three lightest neutrinos. 
The signs in the legend indicate the signs of e za2 , e* a3 , e 1014 = 
±1 for the four possible cases in which CP is conserved. The 
intermediate yellow region is allowed only in the case of CP 
violation. 


Z [eV] 

FIG. 14. Value of the effective Majorana mass \mpp\ as 
a function of sum of the three lightest neutrino masses S in 
Eq. ( |23| ) in the case of 3+1 mixing with Normal Ordering of 
the three lightest neutrinos. The legend is explained in the 
caption of Fig. [12] 
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10“ 2 10" 1 1 
mp [eV] 

FIG. 15. Comparison of the 3cr allowed regions in the mp- 
\mpp\ plane in the cases of 3u and 3+1 mixing with Normal 
Ordering of the three lightest neutrinos. 



E [eV] 


FIG. 16. Comparison of the 3cr allowed regions in the E- 
\mnpp\ plane in the cases of 3v and 3+1 mixing with Normal 
Ordering of the three lightest neutrinos. 


A. Normal Ordering 

Figure [TT] shows a comparison of the best-fit value and 
the lcr, 2 a and 3 cr allowed intervals of the partial con¬ 
tribution /z 4 as a function of the lightest mass m\ with 


those of /i i, (i 2 i M 3 ? which are slightly different from those 
in Fig. [T] because of the contribution of $14 in Eqs. (30)- 
( |32| . One can see that in the 3+1 case it is not pos¬ 
sible to get a total cancellation of \mpp\ in the interval 
mi « (2 — 7) x 10 -3 eV as in the case of 3u mixing (see 
Tab. [I]), because in this interval of mi the contribution 
of M 4 is dominant. However, there is a range of higher 
values of mi between about 0.02 and 0.2 eV in which 
M 4 and Mi have similar values, leading to a possible to¬ 
tal cancellation. For mi > 0.2 eV a total cancellation is 
again not possible because the contribution of Mi is dom¬ 
inant. This behavior is confirmed by Fig. |2j where one 


can see that the value of (see Eq. (|13|)) correspond¬ 


ing to the best-fit values of the partial mass contributions 
is negative for 0.035 < mi <0.1 eV. 

Figure [l2] shows the allowed values of \mpp\ as a func¬ 
tion of mi at different confidence levels. The correspond¬ 
ing intervals of mi for which there can be a total cancel¬ 
lation of \mpp\ are given in Tab. |ll| 

In Fig.[l2]we have plotted separately the allowed bands 
for the eight possible cases in which CP is conserved 
(<U 2 ,o; 3 ,q ;4 = 0 , 7 r), that are the extreme cases which 
determine the minimum and maximum values of | mpp\. 
The areas between the CP-conserving allowed bands cor¬ 
respond to values of \mpp\ which are allowed only in the 
case of CP violation. Unfortunately, these areas are visi¬ 
bly smaller than those in Fig. [3] in the case of 3v mixing. 
This is due to the relatively large uncertainty of M 4 , which 
can be seen clearly in Fig. [II] This uncertainty broadens 
the allowed bands corresponding to the CP-conserving 
cases, leaving little intermediate space. In any case, even 
if the uncertainty of M 4 will be reduced in the future, 
there cannot be a region which is allowed only in the 
case of CP-violation for m\ < 10 -2 eV, where M 4 is dom¬ 
inant and the CP-violating phases are irrelevant. In fact, 
ah the best-fit CP-conserving curves have approximately 
the same value for mi < 10 -2 eV. 


As in the case of 3v mixing, the plot in Fig.[l2]of \mpp\ 
as a function of the lightest mass mi is useful because 
it gives a clear view of the different possibilities for the 
value of \mpp\, but in practice it will be very difficult to 
determine experimentally an allowed region in this plot 
because of the difficulty of measuring the value of the 
lightest mass. Therefore, we calculated also the allowed 
regions in the m^-lm^l and Yr-\mpp\ planes shown in 
Figs. p~3| and [H] w ith the quantities mp and E defined in 
Eqs. \22\ and (|23|) as in the case of 3 v mixing in terms of 
the three standard neutrino masses only. The reason of 
this choice is that mp and E are measurable quantities 
also in the 3+1 scheme. Indeed, considering /3 decay, 
mp quantifies approximately the deviation of the end¬ 
point of the electron spectrum due to neutrino masses 
smaller than the experimental energy resolution [7UH55] . 
whereas the effect of the larger mass m 4 is a kink of the 
Kurie function (see Ref. CD]). In cosmology, the effects 
of the larger mass m 4 can be disentangled from those of 
the smaller masses, because z / 4 becomes non-relativistic 
shortly after matter-radiation equality, much earlier than 
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Vi, u 2 , v 3 . Moreover, it is possible that the contribution 
of 7714 to the energy density of the Universe is suppressed, 
for example by a large lepton asymmetry or an 

enhanced background potential due to new interactions 
in the sterile sector [Ml [95] , or a larger cosmic expansion 
rate at the time of sterile neutrino production [M], or 
MeV dark matter annihilation [97] . 

From Figs. p~3] and 14 one can see that the intervals of 
171/3 and E for which there can be a complete cancellation 
of the three partial mass contributions to \mp p\ (given 
in Tab. [II| are much larger than those in Figs. [4] and [ 5 ] 
for the standard Zu mixing case, and \mpp\ > 0.01 eV for 
any value of the unknown phases 0:2, 0:3, ck 4 only for the 
relatively large values mp > 0.25 eV and E > 0.8 eV. 

It is useful to compare the allowed regions mp-\mpp\ 
and T,-\mpp\ planes obtained in the cases of Zv and 3+1 
mixing with Normal Ordering of the three lightest neu¬ 
trinos. Figures [15] and [16] show this comparison for the 
3 a allowed regions. One can see that, if the Normal Or¬ 
dering will be established by oscillation experiments (see 
Refs. DEI), with measurements of mp and \mpp\ and/or 
E and \mpp\ it may be possible to distinguish Zv mixing 
and 3+1 mixing if the measured values select a region 
which is allowed only in one of the two cases. It is in¬ 
teresting that there are two regions allowed only to 3+1 
mixing: one with \mpp\ smaller than that in the case of 
Zv mixing and one with \mpp\ larger than that in the case 
of Zv mixing. At least a part of the second region is acces¬ 
sible to the next generation of neutrinoless double-beta 
decay experiments (see Refs. [36b4T] ). The only /3-decay 
experiment under preparation with the aim of exploring 
the sub-eV region of mp is KATRIN [98], which will have 
a sensitivity of about 0.2 eV that is not sufficient to ex¬ 
plore the upper part of the region in Fig. [15] allowed only 
in the case of 3+1 mixing. On the other hand, cosmo¬ 
logical observation may be able to measure the sum of 
the three light neutrino masses down to the lower limit 
of about 5.6 x 10 -2 eV l99l. 


B. Inverted Ordering 

The best-fit value and the lcr, 2cr and 3cr allowed inter¬ 
vals of the partial mass contributions to \mpp\ in the case 
of 3+1 mixing with Inverted Ordering of the three light¬ 
est neutrinos are shown in Fig. [XT] One can see that there 


TABLE III. Ranges of m 3 , mp and E for which there can be 
a complete cancellation of the four partial mass contributions 
to \mpp\ for the best-fit values (b.f.) of the oscillation param¬ 
eters and at lcr, 2cr and 3cr in the case of 3+1 mixing with 
Inverted Ordering of the three lightest neutrinos. 

b.f. lcr 2cr 3cr 

m 3 [10~ 2 eV] < 9.1 < 11.4 < 15.5 < 19.3 

mp [1(T 2 eV] 4.8 - 10.3 4.8 - 14.2 4.7 - 17.6 4.6 - 22.5 
E [HT 2 eV] 9.8 - 29.9 9.7 - 41.8 9.5 - 52.3 9.4 - 67.1 



10~ 4 10“ 3 10 “ 2 10 _1 1 


Lightest mass: m 3 [eV] 


FIG. 17. Best-fit values (b.f.) and lcr, 2cr and 3cr allowed 
intervals of the four partial mass contributions to \mpp\ in 
Eq. (29) as functions of the lightest mass m 3 in the case of 3+1 
mixing with Inverted Ordering of the three lightest neutrinos. 


1 


10” 1 


> 

Cl) 



ca 

CCL 

E 


1 



Lightest mass: m 3 [eV] 

FIG. 18. Value of the effective Majorana mass \mpp\ as a 
function of the lightest neutrino mass m 3 in the case of 3+1 
mixing with Inverted Ordering of the three lightest neutrinos. 
The legend is explained in the caption of Fig. [12] 


can be a total cancellation between the partial contribu¬ 
tion 114 and the dominant in the case of Zv mixing (see 
Fig. [ 6 ]) for m 3 < 0.1 eV. Indeed, Fig. [ 8 ] shows that the 
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FIG. 19. Value of the effective Majorana mass \mpp\ as 
a function of sum of the three lightest neutrino masses E in 
Eq. ( |23| ) in the case of 3+1 mixing with Inverted Ordering of 
the three lightest neutrinos. The legend is explained in the 
caption of Fig. [12] 



FIG. 21. Comparison of the 3cr allowed regions in the mp- 
1 771 / 3 / 3 1 plane in the cases of 3v and 3+1 mixing with Inverted 
Ordering of the three lightest neutrinos. 



E [eV] 


E [eV] 


FIG. 20. Value of the effective Majorana mass \mpp\ as a 
function of effective electron neutrino mass mp in Eq. (22) in 


the case of 3+1 mixing with Inverted Ordering of the three 
lightest neutrinos. The legend is explained in the caption of 
Fig. pl 


value of mfpp (see Eq. jljj )) corresponding to the best-fit 


FIG. 22. Comparison of the 3cr allowed regions in the E- 
1777/3/31 plane in the cases of 3v and 3+1 mixing with Inverted 
Ordering of the three lightest neutrinos. 


values of the partial mass contributions is negative for 
7773 < 0.09 eV. 

Figure [18] depicts the allowed regions in in the 7773 - 
1777/3/31 plane. Comparing Fig. |l8| with Fig. [7|one can see 
that the predictions for \mpp\ are completely different in 
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the Zv and 3+1 cases if there is an Inverted Ordering 
of the three lightest neutrinos, in agreement with the 
discussions in Refs. [ 13] H5IH3] . The ranges of values 
of m 3 for which there can be a complete cancellation of 



in the cases of Zv and 3+1 mixing with Inverted Order¬ 
ing of the three lightest neutrinos. If the Inverted Or¬ 
dering will be established by oscillation experiments (see 
Refs. PH]), it will be possible to exclude Zv mixing in 
favor of 3+1 by restricting mp and \mpp\ or X and \mpp\ 
in the corresponding large region at small \mpp\ allowed 
only in the 3+1 case. 


IV. CONCLUSIONS 

We have presented accurate calculations of the effec¬ 
tive Majorana mass \mpp | in neutrinoless double-/^ decay 
in the standard case of Zv mixing and in the case of 3+1 
neutrino mixing indicated by the reactor, Gallium and 
LSND anomalies (see Refs. [HHH]). We have taken into 
account the uncertainties of the standard Zv mixing pa¬ 
rameters obtained in the global fit of solar, atmospheric 
and long-baseline reactor and accelerator neutrino oscil¬ 
lation data presented in Ref. [3] and the uncertainties 
on the additional mixing parameters in the 3+1 case ob¬ 
tained from an update [34] [35] of the global fit of short- 


baseline neutrino oscillation data presented in Ref. [RJ] . 

We have shown that the predictions for \mpp\ in the 
cases of Zv and 3+1 mixing are quite different, in agree¬ 
ment with the previous discussions in Refs. [13] I25b33] . 
Our paper improves these discussions by taking into ac¬ 
count the uncertainties of all the mixing parameters and 
presenting all the results at la, 2cr and 3cr. 

We have presented accurate comparisons of the allowed 
regions in the planes 777 ,^-lm^l and E-lm^l of measur¬ 
able quantities, taking into account the two possibilities 
of Normal and Inverted Ordering of the three light light¬ 
est neutrinos. We have shown that future measurements 
of these quantities may distinguish the Zv and 3+1 cases 
if the mass ordering is determined by oscillation experi¬ 
ments (see Refs. PH]). 

We have also introduced in Section |II A| a relatively 
simple method to determine the minimum value of \mpp\ 
in the general case of N -neutrino mixing. 
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